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ABSTRACT 

Context. We present the results of new spectral diagnostic investigations applied to high-resolution long-slit spectra 
obtained with the Hubble Space Telescope Imaging Spectrograph (HST/STIS) of the jet from the T Tauri star RW Aur. 
Aims. Our primary goal is to determine basic physical parameters (electron density no and electron temperature Te, 
hydrogen ionisation fraction Xc, total hydrogen density riH, radial velocity and the mass outflow rate Afj) along both 
the red- and blueshifted lobes of the RW Aur jet. 

Methods. The input dataset consists of seven long-slit spectra, of O'.'l spatial resolution, taken with the STIS slit parallel 
to the jet, and stepped across it. We use the Bacciotti & Eisloffel (1999) method to analyse the forbidden doublets 
[OI] AA6300,6363, [SII] AA 6716,6731, and [Nil] AA6548,6583 A to extract rie, Te, Xe, and nn. 

Results. We were able to extract the parameters as far as 3'.'9 in the red- and 2'.'l in the blueshifted beam. The electron 
density at the base of both lobes is close to the critical density for [S II] emission but then it decreases gradually with 
distance from the source. The range of electron temperatures derived for this jet (To=10^-2 x lO** K) is similar to those 
generally found in other outflows from young stars. The ionisation fraction Xc varies between 0.04 and 0.4, increasing 
within the flrst few arcseconds and then decreasing in both lobes. The total hydrogen density, derived as nn ~ Uc/xc, 
is on average 3.2 x lO^cm"^ and shows a gradual decrease along the beam. Variations of the above quantities along the 
jet lobes appear to be correlated with the position of knots. Combining the derived parameters with Vy measured from 
the HST spectra and other characteristics available for this jet, we estimate Mj following two different procedures. The 
mass-outflow rate Mj is moderate and similar in the two lobes, despite the fact that the well-known asymmetry in the 
radial velocity persists close to the source. Using the results of the BE diagnostics we find averages along the first 2'.'l 
of both fiows (a region presumably not yet affected by interaction with the jet environment) of 2.6 x 10~^ Mq yr~^ for 
the red lobe and 2.0 x 10~^ Mq yr~^ for the blueshifted fiow, with an uncertainty of ±log A/q = 1.6. 
Conclusions. The fact that the derived mass outfiow rate is similar in the two lobes appears to indicate that the central 
engine is constrained on the two sides of the system and that the observed asymmetries are due to environmental 
conditions. Possible suggestions for the origin of the differences are discussed. The RW Aur jet appears to be the second 
densest outfiow from a T Tauri star studied so far, but its other properties are quite similar to those found in other 
jets from young stars, suggesting that a common acceleration mechanism operates in these sources. 

Key words. ISM: jets and outfiows - Stars: circumstellar matter - Stars: pre-main sequence: - Stars: mass-loss - Stars: 
individual: RW Aur 



1. Introduction 



Highly-coUimated jets from yo ung stellar objects (YS Os) have been thought for many years to be an essential ingredient 
of the star for mation process (|Mundt fc Friedlll983l) . Their origin, ejection, and coUimation is, however, still a matter 
of debate (e.g. iRav et al.l[2007l ). Today, high-resolution spectroscopy from space and the ground allows us new insights 
into the physical conditions of these jets close to their source. This opens up new ways to compare observations with 
the different magncto-hy drodynam ic (MH D) models propose d for the jet launch and acceleration (e.g. ICamenzin dl ll990l : 
iKonigl fc Pud ritz 2000: Shu et all , 2000: Matt et al.l l2003: F erreira fc Cassell2004[) . During the past few y ears, high reso- 
lution studies of several jets fr om young s tars on subarcsecond s cales have become available - e.g HH 30 ( Bacciotti et al.l 



lution stuaies ot several lets tr om young s tars on subarcsecona s cales nave become available - e.g tin 60 
119991: iHartigan fc Morsell2007l). DG Tau (iBacciotti et al.|[2000l: iLavallev-Fo uguet et all 120001: iPvo et al 



20031) . HL Tau 



(|Pvo et al.ll2006f ). RW Aur (jPvo et al1l2006D . and LkHa 233 (jPerrin fc Grah am 20071: iMelnikov et al.ll2008f) . The present 



study is a continuation of the spectroscopic analysis of RW Aur p erformed with the S pace Telescope Imaging Spectr ograph 
(STIS ) on-board the Hubble Space Telescope (HST), started bv lWoitas et al.l (|2002l hereafter Paper I), and iWoitas et al.l 
(l2005l) . 



* Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. 
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The classical T Tauri star RW Aur is a double system in the Taurus- Auriga star forming region. The primary. 
RW Aur A (Sp. K1-K4), is one of the optically brightest T Tauri stars with V = 10^1. The secondary, RW Aur B 
{V = 12™7), is located at a projected separation of 1'.'2 and a position angle of 256° with respect to the primary. 
RW Au r might be a mult iple hierarchical system, although this is not clear. RW Aur A is a s uspected spectros copic 
binary (jPetrov et al.ll200lt ). Another component , which was resolved very close to the secondary (jGhez et al.lll993[) . was 
probably a false detection (jWhite fc Ghej|200l[ ). 

A bipolar jet, with well-coUimated blue- and redshifted lobes emanating from RW Aur A (HH229), was found by 
iHirth et al.l ()1994f ). Further w ork traced the emission out to 15" in the redshifted lobe (jPougados et al.l 120001) and over 
100" in the blueshifted beam (jMundt fc Eisloffelll998t ). Radial velocity measurements of several forbidden lines revealed 
an interesting asymmetry: the radial veloci ty of the blueshift ed lobe reaches about — 190kms~^, while a velocity of only 
+ 100kms~^ is measured on the red side (jHirth et al.l 1199-1 ). Another unusual prope rty of this jet is that fo r the first 
10" the redshifted jet is brighter (in the red [SII] doublet) than the blueshifted beam (|Mundt fc Eisloffeilll998[ ). whereas 
in general blueshifted jets are brighter. Beyond this distance, the redshifted RW Aur jet strongly drops in brightness, 
becoming much fainter than the blueshifted flow. 

The physical parameters along the RW Aur jet were first studied bv lBacciotti et all (Il996f), using sp ectra taken from 
the ground at moderate spatial and spectral resolution, and subseguentlv bv iDougados "et all (|2000f ). again from the 
ground but at higher spatial resolution using adaptive optics. These estimates are discussed in detail later in the present 
paper. A substantial advance in our understanding of this jet came, however, with the observations at high angular and 
spectral resolution performed with HST/STIS. In this case, seven long-slit spectra of the RW Aur jet were ta ken at O'.'l 
spatial resolution, with the STIS slit parallel to the jet, and stepped every 0"7 across it. In lWoitas et al.l ()2002f ) (hereafter 
Paper I) we presented reconstructed velocity channel maps of the bipolar jet in the forbidden lines, giving indications 
of the morphology and kinematics of the jet close to the star and at high resolution. In the second paper of the series, 
IWoitas et al.l (|2005f ). we reported signatures of rotation in both jet lobes within the first 1"5 from the central source 
(about 200 AU at the distance of RW Aur). Similar sign atures of rotation have been discovered in other jets from T 
Tauri stars (jBacc iotti et al."2002': 'Coffev et al."2004'. '2007^. Such observations support the magneto-centrifugal launching 
scenario for jets (cf Ferreira et al. 2006. : .Pudritz et al. 2007t, and references therein), which predicts a transport of angular 
momentum from the disk to the jet by the accretion/ejection 'eng ine'. Recently, how ever, the disk surrounding RW Aur 
has been found to rotate in the opposite sense of the bipolar jet (jCabrit et al.ll20d6[ ). An explanation for this apparent 
discrepancy has not yet been found. 

In this third paper of the series, we continue the analysis of our high- resolution HST/STIS spectra of the RW Aur 
jet, to derive the physical quantities in the emitting gas from the observed line ratios. Applying specialised diagnostic 
techniques we estimate the electron density and temperature, the ionisation fraction, the total hydrogen density and 
the mass outflow rate along the jet. Here we report on such results in detail, and we offer a new analysis of the spatio- 
kinematical structure of the bipolar jet extracted from the high-resolution position-velocity maps. 

The paper is structured as follows: the observations and details of data reduction and analysis are described in Sect.j^l 
In Scction[3]wc illustrate the estimates of the physical quantities along both lobes of the jet. Then in Section[4]we contrast 
the derived jet properties of the opposing lobes, and we compare the estimates with values derived previously and in 
other T Tauri jets. Finally, we summarize our conclusions in Sect. |SJ 



2. Observations and Data Analysis 

Our long-slit HST/STIS spectra of the RW Aur jet were ta ken on 10 Dec 2000. Obser vational details and the basic 
steps of the data reduction have been described in iPaperl and in IWoitas et al.l (|2005[ ). Briefly, a slit of width O'.'l, 



length 52"and parallel to the RW Aur jet axis (P.A. = 130°. lDougados et al.ll2000l i was moved across the jet in steps of 
0'.'07 from south-west (slit SI) to north-east (slit S7). Thus, our dataset consists of seven spectra (S1,..,S7) with a spatial 
sampling O'.'05/pixel, and, with the chosen grating (G750M), a spectral sampling of 0.554 A/pixel, corresponding to about 
25kms~^/pix in the selected wavelength range (6300 - 6860 A). The basic reduction of the spectra has carried out by the 
HST/STIS pipeline, and includes flatficlding, cleaning of cosmics and bad pixels, and calibration in surface brightness and 
wavelength. The spectra include prominent forbidden hues - [01] AA6300,6363, [Nil] AA6548,6583, [SII] AA6716,6731, in 
addition to the stellar continuum and strong Ha emission. 

In order to trace the jet in the forbidden lines as close as possible to the source, the continuum has to be removed. 
This is easily done by subtracting the continuum fitted row by row with a high-order polynomial on a sampling region 
that excludes the emission lines. Ha is also subtracted in the same way, as it is not usable for the diagnostics and its 
high luminosity affects the [Nil] lines. Subtraction residuals, however, grow rapidly for positions close to the source, so 
the area inside ±0'.'15 is excluded from further analysis and it is masked in all the figures presented in this paper. 

In order to improve the quality of the extraction of the flux in the lines and to calculate statistical errors, we determined 
the background from running averages of boxes of size 7x8 (170kms^^x 0'.'4) pixels moving along both sides of the 
lines in each individual slit, and we then subtracted it from the fluxes of the emission lines along the jet lobes. 

An example of the cleaned spectra in each of the seven slit positions for the [S II] A6731 line is provided in Fig. [TJ 
in the form of Position- Velocity (PV) diagrams. Similarly, PV diagrams for each emission line, this time obtained from 
coadding the spectra obtained in all the seven slit positions, are shown in Fig. |21 The PV maps approximately span 
velocities in th e range from —300 to -f 250kms~-'^, and are corrected for the heliocentric radial velocity of RW Aur A 
f-l-23.5kms~^. IWoitas et"ani2005t ) . The contours are all scaled to the same range of intensity in each figure. 
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To apply the diagnostic technique, the resulting pure forbidden flux was further integrated over the line width. This 
step was dictated by the fact that the emission is concentrated in a highly coUimated jet moving in a narrow range of 
radial velocities (see iPaper J ). Therefore, the signal is not strong enough at the lateral borders of the jet, and at the 
extremes of the velocity range, to allow a meaningful application of the diagnostic technique at the original spatial and 
spectral resolution. High angular resolution (O'.'l on two pixels) is however retained along the jet. 

To perform the diagnostic analysis, line ratios were then calculated every 0'.'05 along the jet axis and used to derive 
the physical quantities of interest (electron density Hq, hydrogen ionisation fraction Xc, el ectron temperature Te, and 
total hydrogen density tih along the jet). To this aim, we followed the method described in lBacciotti fc Eisloffell (|1999( ) 
(hereafter BE99). The so-called BE technique illustrated in that paper requires as input the flux ratios [SII] A6716/A6731, 
[01] A(6300+6363)/[SII] A(6716+6363), and [Nil] A(6548+6583)/[OI] A(6300+6363) (in the following the latter two ra- 
tios will be referred to simply as [01] /[SII] and [Nil]/ [01]). The comparison between the ratios calculated theoretically 
and the ones observed gives es timates of the physical quantities. The ratio of the [S II] doublet is used for the determi- 
nation of ric (|Osterbrocklll989( ). up to densities of about 2.5 x lO^'cm"'^, above which the ratio reaches a plateau. The 
ratio [Nil]/ [01] mostly depends of the ionisation fraction .Tq, and [01] /[SIX] is dependent on both Xc and Te. 

The errors on the physical quantities extracted with the BE technique originate from the background noise, (calculated 
from the RMS of the fluxes in running average background boxes) but also come from other sources. Whereas the error of 
the electron density ric depends mainly on the error of the ratio [SII] A6716/A6731 propagated from the noise values, the 
estimates of temperature Te and hydrogen ionisation fraction Xo depend both on the variations of the ratios [01] /[SII] 
and [N II] / [O I] respectively and on the uncertainty of the Ug estimates. Another source of error is the uncertain extinction 
Av not only at the source, but also along the jet. Nevertheless the diagnostic analysis is based on ratios of lines quite close 
in wavelength, and thus the results are almost independent on dereddening correction, at least for extinction typical of T 
Tauri stars. In order to evaluate the impact of extin ction in the RW Aur c ase, however, we used a plausible range of 
values from (P03 to CP72, derived from the study of lWhite fc Ghed (|200l[ ). to correct forbidden Hue fluxes. Comparison 
of the resulting difference with the noise values shows that the former are only a few percent of the S/N errors. In turn, 
we expect that errors due to variation of extinction along the jet beams are negligible, and therefore we used a constant 
average value A^ = 0™4 everywhere along the jet to find the physical parameters. As we will discuss in Section !?^ 
however, extinction variations may play an important role in calculations involving the absolute flux of the lines, as 
used to estimate the jet mass outflo w rate. The results of the BE technique also depend on the adopted abundance of 
chemical speci es. Podio et al.l (j2006[ ) showed that the interstellar gas composition determined for the Orion nebula by 
lEsteban et al.l (|2004f) led to the most reliable diagnostic results for jets in various regions. Thus, we adopted the same 
abundance set, assuming RW Aur star formation region has been exposed to the same enrichment processes. 

Finally, we note that a fundamental assumption for the application of the BE technique is that the minimum region 
sampled is homogenous in the physical conditions of the plasma, in such a way that the excitation of the different lines 
can b e well described by a unique electron temperature, electron density and ionisation fraction. As demonstrated in 
lBE99l this condition holds in the cooling region of shock fronts with shock velocities between 30 and 70kms~-'^, typical 
of jets. In the present case, however, the line ratios are calculated after coadding the spectra in the transverse direction, 
which implies the additional assumption that at each distance from the star the gas conditions are nearly homogenous in 
the jet transverse direction. In the present case, such an assumption is justifled by the fact that the images and spectra 
of the flow show an extremely well collimated jet with very little dispersion in velocity. This supports a scenario in which 
the jet flow intersects the internal bow shock fronts in a small area at their apex, in which the front can be thought 
to be perpendicular to the flow and cover the jet transverse section with similar shock strength. In turn, this indicates 
homogenous excitation conditions of the gas across the transverse jet section. 

In order to calculate the mass outflow rate in the jet one needs an accurate determination of the jet velocity. Spatio- 
kinematic information (radial velocity maps including peak radial velocities) was derived from the spectra every Qf.'b from 
the star with Gaussian fits after correction for uneven slit illumination (see Section [ 



3. Results 

3.1. Preliminary inspection of the line ratios 

The distribution of the line ratios together with their errors allows us to estimate the data quality of the input values. 
This, in turn, gives an indication of the reliability of the derived physical parameters at each position from the source 
and suggests ways to improve the calculation. 

In the top panels of Figs. [3] and 2] we report the values of the line ratios [SII] A6716/A6731, [01] A6300/A6363 and 
[Nil] A6583/A6548 in the red and blue lobe, respectively, compared with the global intensity profile along the jet. The 
ratios and the profile are calculated from the lines extracted from the original spectra with a spatial sampling of O'.'OS, 
and subsequently integrated in radial velocity and across the flow. In the plot for the [S II] ratio the limiting values for 
the He diagnostics of 0.465 (high density limit, corresponding to the critical density for coUisional de-excitation) and 
1.445 (low density limit) are shown. Values outside this range do not provide a good diagnostic determination for the 
electron density, and in the corresponding positions we have assumed that the electron density is equal to the critical 
density, 2.5x 10^ cm~'^, or 50cm~^, respectively, thus giving a lower or upper limit to the actual electron density (see 
next Section). 

In addition, the [O I] A6300 / A6363 and [N II] A6583/ A6548 ratios should be equal to the ratio of Einstein's A coefficients 
for spontaneous emission in each species, which is constant and equal to about 3 in both cases. In the positions where 
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values of these ratios show consistent deviations from the theoretical values, either one or both lines of the doublet were 
weak and/or corrupted by imperfect subtraction of the stellar continuum or isolated hot pixels, or by the presence of 
strong Ha emission. The situation is more critical in both lobes for the [Nil] A6583/A6548 ratio: in the first 0'.'5 from 
the source this ratio is often exceedingly low, probably because of an underestimate of the [N II] A6583 line flux caused 
by oversubtraction of the red wing of the Ha line, that is very strong close to the source. Farther from the star, and 
especially after 2", the situation is further complicated by the faintness of the [Nil] A6548 line, that, after subtraction of 
the Ha halo, results in recognizable regions of artificial brightening or depression of the line, leading to uncertain and, in 
some cases, even negative values of the fiux. The negative values are clearly unphysical and have been either dropped or 
corrected setting [Nil] A6548 = 1/3 [Nil] A6583 in the corresponding positions (identified by filled symbols in Figs. [3] and 
131). In the other locations where the departures are observed it is difficult to establish which of the two lines is corrupted, 
so we preferred to keep the measured values, using the panels of Figs.[3]and|3]as a flag map to defective positions. 

We note also that in the blue lobe, the blueshiftcd [01] A6300 line lies close to the edge of the CCD (see Fig.[2|), and 
the emission on the left border is partly cut. Nevertheless, by comparing at each distance from the star the integrated 
flux of the [O I] A6300 spectral image with the integrated flux from a similar image obtained from multiplying by 3 the 
[01] A6363 (Figs. [5]), we have checked carefully that the flux lost is negligible everywhere (not more than a 0.5% of the 
integrated flux), which is also confirmed by the plot of the [01] A6300/A6363 ratio. 

Finally, in the bottom panels of Figs. [3] and d] we show the ratios [Nil]/ [01] and [OI]/[SII] that are used as input 
to the diagnostic code. In both lobes the [Nil]/ [01] ratio presents a slight monotonic increase with distance from the 
source, while the [OI]/[SII] ratio shows a pronounc ed decrease from the positions close to the source toward the je t 
beam. The latter effect, already reported for HH30 (jBacciotti et al.l [19991 ) and other HH flows (jHartigan et al.l [19951 ). 
occurs in the high density regions close to the jet source, where the [S II] lines are quenched, while the [O I] lines are not 
affected, having higher critical density (~ 2 x 10^ cm^"^). In both ratios, local enhancements are seen in the red lobe at 
the position of knots J3 and J4. 

The scatter with respect to the general trend is larger in locations corresponding to minima of the emission flux, as. 
e.g., in the regions between the knots (e.g. at 1'.'3 and 2'.'8 in the red lobe). Spatial binning over the original sampling 
of O'.'05/pixel would have increased the signal-to-noise in these portions, but we avoided its use in order not to loose the 
unique information coming from the high angular resolution of HST data. 

Another critical region is the first 0'.'2 close to the star, because of imperfections in the subtraction of the continuum: 
despite the line ratios have been analysed here, this zone is excluded from the subsequent diagnostic analysis. Similarly, 
in both lobes the faint final portions of the beams (after 3'.'8 and 2" in the red and blue lobes, respectively) could not be 
adequately treated by the BE technique and are not reported in Figures^ HI and subsequent figures. 

3.2. Physical quantities along the flow 

The set of line ratios illustrated above constitutes the input for the spectral diagnostics that leads to the determination 
of the physical quantities in the emittin g gas. The calcu lation is performed numerically using a code developed to apply 
the BE technique to large datasets (see iMelnikov et aI1[2008il . The physical parameters derived for the two lobes of the 
RW Aur jet are presented in graphical form in Figures [6] and [71 The main features are described below, separately for 
each lobe. The detailed kinematics of the flow is then illustrated in Sect. 13.31 Globally, the implications of these results 
are discussed in Sect. SI 

3.2.1. Red lobe 

In Fig. [6] the top panel shows, for reference, the distribution of the coadded intensity in the forbidden lines integrated 
in velocity and across the beam, along the first 4" from the source with a spat i al sam pling of O'.'05/pixel. The identified 
knots are labelled J1-J6, follo wing the naming sch e me of [Mundt fc Eisl5ffel (|1998[ ) (the three brightest knots at the 
same jet section are labelled bv iLopez-Martm et al.l (|2003l ) as R4-6). Similarly, we report the profile of the peak radial 
velocity in [SII] A6731 and [01] A6300 for reference. This quantity is illustrated in more detail in Fig. [H] and analysed in 
Sect. [gX^ 

The next panel presents the electron density, rie with its propagated error from the [S II] lines. In the first 0'.'5 from 
the star Uc reaches the highest values, just below the critical density for [SII] diagnostics, richer = 2.5 x 10'*cm~'^ (error 
bars without the upper/lower tick indicate that errors on the [SII] lines lead the ratio above the limiting value). Then 
He gradually decreases, with an average value of 3400 cm^"^ along the beam. A local minimum is found at 1'.'4 and 
pronounced scatter is visible between 2'.'5 and 3'.'1, coincident with minima of the line emission in the top panel. 

The hydrogen ionisation fraction Xa (fourth panel) gradually increases along the flow, from about 0.02 to an average 
of 0.07, but it reaches a value as high as 0.2 at 1'.'5, on the left edge of the ,13 knot. Other local maxima are present at 
the position of knots J2, J4, and J6, while the scatter in the faint portion between J5 and J6 is probably due to the 
uncertainty introduced by weak [Nil] emission. It is interesting that at knot J3 the maxi mum of Xg is actua l ly at the 
upstream edge of the structure. Such behaviour has also been noted in the HH30 knots bv lHartigan fc Morsel (|2007l ). 

The electron temperature To calculated from the BE technique is illustrated in the next panel. The base of the jet 
presents the highest temperature, at ~ 1.5 x 10'*K. This gradually decreases to 9000 K at 1'.'3, but then it sharply rises 
back to an average of ~ 1.5 x 10* K at the position of knots J3 and J5 (between 1'.'5 and 2'.'2). In the following faint 
section of the jet the average is quite moderate, ~ 10** K, and progressively decays toward the end. 
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The second to last panel from top illustrates the total hydrogen density nn, derived as 71^1 x^. The errors propagated 
from the other quantities are large in many positions, but on average one can recognize a slow decrease with distance 
from the source from log(nH) ~ 5.7 close to the source to ~ 4 at the last position examined. In the region of knots J3-J5, 
however, the density appears to locally follow the opposite trend, and likewise at knot J6. 

The bottom panel of Fig. [H] shows the mass outflow rate. The derivation will be discussed in detail in Sect. 13.41 and 
using Fig. [TOl 

3.2.2. Blue lobe 

In this section we analyse the results for the blueshifted lobe, which is fainter than the redshifted one over the examined 
distances from the source (2"). The physical quantities derived in this region with the BE diagnostics are presented in 
Fig.H 

The distribution of the radiative flux along the beam in the examined lines is shown in the top panel. In contrast to 
the red beam, the blue lobe does not present a prominent maximum close to the star, but it drops steadily from the first 
examined positions. The distribution of forbidden inte nsity is flatter, wi t h less prominent knot feat ures at 0^.^5 and 1"2, 
marked as All and A12 following the nomenclature of iMundt fc Eisloffell (|l998( ) (or B3, B4 foUowing lLopez-Martfn et al.l 
|2003[ ). The radial velocity in the second panel from top is nearly constant at — 180kms~^ up to knot A12, and then it 
decreases to — 160kms~^ at 2'.'0. Note, how velocities are much higher, in absolute terms, than the corresponding ones in 
the red lobe. The kinematic properties are presented in more detail in Sect. 13.31 while the velocity asymmetry is discussed 
in Sect.H 

At the base of the lobe (first 0'.'4 from the source), the electron density (third panel from top) is higher than in the 
red lobe, with ?io > rtcr in various positions (marked by arrows). Beyond this distance, the trend of the electron density 
is not uniform: rto drops gradually, but there is also a local maximum at 1'.'2, and a slight rise of tIo at the end of the 
section (at about I'.'S). Within the examined 2'.'1 the average value of log(nc) is 3.56, which is almost the same as in the 
red lobe (log(ne) = 3.61) within similar distances from the source. 

The ionisation fraction (fourth panel) shows a substantial scatter along the way to 2'.'1. Nevertheless, there is 
marginal evidence that Xf. increases with distance from the source reaching local maxim a of 0.3 and 0.4 at the positions 
of th e knots. Such en hancement of x^ is also observed in the fainter lobes of other jets ijBacciotti et al.l [19991 : iBacciottH 
I2002t iHartigan fc Morse 2007 ). Beyond 1'.'2, Xe appears to drop to 0.1, but the errors are higher here. The average value 
Xe for the whole section is 0.23, considerably higher than in the redshifted jet [x^ = 0.08). 

The electron temperature To (fifth panel) is also higher than in the red lobe, being on average To = 1-6 x 10"* K over 
the studied section. Despite the large dispersion of the values, a slight increase can be identified in the first 0'.'6; followed 
by a local minimum at 10"* K. After that, Tq varies around 1.6 x 10^ K which coincides, within the errors, with the mean 
value. 

In the following panel one finds the total hydrogen density nn, calculated as nic/xc- This shows a gradual decrease 
along the beam, with log(nH) varying from 5.0 to 3.6. The arrows refer to positions in which Uc is above the critical 
density, and thus represent lower limits. In general, there are no large deviations from the mean value of log(nH) = 4.4. 

Finally, in the bottom panel of Fig. [7| we report our determination of the mass outflow rate. Again, its derivation will 
be discussed in Sect. 13.41 

3.3. Kinematic properties of the flow 

The seven HST/STIS spectra of the RW Aur jet, taken at high spatial and spectral resolution, provide the most detailed 
description of the kinematics of a bip olar coUi mate d stellar jet c l ose to its so urce to date. This information has been 
described in part in previous papers (jPaper j . and IWoitas et al.ll2003l [2005^ . In this Section we complete the study 
providing an alternative description of the radial velocities along and across the jet lobes. 

General indications come already from the position- velocity (PV) diagrams illustrated in Fig. [21 We recall that these 
isophotic maps, spanning the range from —300 to +250 km s^^, originate from coadding the spectra in the seven slit 
positions, and hence loosing the spatial resolution across the jet. In these maps the radial velocity asymmetry is evident, 
as the material moves at about lOOkms"^ and — 190kms~^ in the red and blue lobe, respectively. The velocity, however, 
presents only little dispersion and variation along each lobe. 

The knots do not present any identifiable double velocity component, but the velocity dispersion varies in the different 
lines. For example, the measured Full Width Half Maxima {FWHMs) show that the jet has the lowest dispersion in [Nil] 
and in the red lobe. Also, the dispersion is larger within knots than between them. On the other hand, the [S II] emission 
appears to be more uniform than the [O I] emission, while [Nil] is mostly concentrated at the knots. It is interesting that 
unlike the other forbidden lines, the flux in the [Nil] doublet appears to be stronger in the blue lobe, testifying to its 
higher excitation (see Sect. 13. 2T^ . 

3.3.1. Profile of the radial velocity along the flow 

In Fig. [5] we illustrate the variation of the peak radial velocity along the jet obtained from Gaussian fits to each emitted 
line separately, integrated across the jet width (that is, from the coadded seven spectra). The integrated lines have 
prominent profiles along most of the beam, which gives a good accuracy to the fit, (uncertainty is about ±5kms~^, 
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excluding the fit to [O I] A6300 (see below) , and except for a number of positions in which the [N II] lines are faint) . The 
first two panels from top show a comparison of the radial velocities in the two lobes on the same scale, while the bottom 
panel illustrates the results for the entire red lobe. In general, the relative variation of the velocity is quite limited in 
both lobes, being less than 15-20% of the average flow velocity all along the flow, and reaching its maximum value in 
the positions corresponding to knots J3-J4. 

In the red beam (top and bottom panels) the trends in the radial velocity of the four bright lines [O I] AA6300,6363 
and [SII] AA6716,6731 are similar, with little variation around the average value of +100kms~^ over the whole section. 
At the beginning of the jet the [01] lines appear to have higher and more constant radial velocity than the [SII] lines, 
whose velocity, however, gradually increases to reach that of the [01] lines at about < 0'.'6 from the source. A possible 
interpretation for this difference may come from their higher critical density. Close to the source the [01] lines trace 
preferentially denser material, probably located closer to the axis, while [S II] emission includes less dense material located 
towards the outer edges of the flow. This would suggest that the gas closer to the jet axis moves faster and may have been 
accelerated more than the gas at the edges. As mentioned above, the radial velocity correlates positively with intensity 
at knots J3-J4, being higher (about 120kms~^) at these locations. This would be in line with the knots representing 
bow-shaped working surfaces in the flow. No clear correlation is observed, however, for the other knots. Despite the [Nil] 
lines being fainter in our HST/STIS spectra, velocity determinations were successful at a number of positions. In general 
the variations of the [N II] radial velocities follow those of the other lines, and reach slightly higher velocit ies at the bright 
knots , consistent with [Nil] emission tracing faster and highly excited material close to the axis (cfr. iBacciotti et al.l 
[2OOOI) . 

In the second panel of Fig. [S] we report the values of the radial velocity measured in the examined portion of the blue 
lobe. One notes again the asymmetry in the radial velocity of the two lobes, present in all lines, that makes the blue lobe 
faster than the red one by 1.5-2 times. This will be discussed in Sect. H) Secondly, in the first arcsecond the [01] A6300 
line presents a small shift with respect to the other forbidden lines, including [O I] A6363. The shift is probably due to 
the fact that [01] A6300 lies close to the CCD border, which causes a lack of accurateness in the Gaussian fit. Values 
beyond 0'.'7, however, are reliable and still show a slightly higher velocity in the [01] lines (as well as the [Nil] lines) 
with respect to the [SII] lines, again probably related to the fact that [01] and [Nil] lines trace denser material closer to 
the axis. The gradual drop after 1" from the source correlates with the decrease in intensity. On the other hand, there 
is no correlation with the position of the knots. This behaviour is quite unusual and will be examined in the Discussion 
section. 



3.3.2. Spatial distribution of radial velocity. 

We have presented above the peak velocity profiles obtained from lines coadded from the seven spectra. From our dataset, 
however, one can resolve spatially the velocity field in two dimensions (along and across the jet). To this aim, spatial 
maps of the peak radial velocities have been constructed. First, to increase the S/N ratio we coadded in each spectrum 
the two components in each of the [01], [Nil], and [SII] doublets, with the exception of [01] A6363 in the blue lobe, as 
the other oxygen line lies on the border of the CCD and its Gaussian fit proved to be critical. Then for each spectrum, i.e. 
each slit position, and at each distance from the star we performed a Gaussian fit of the spectral profiles and found the 
peak velocity. We adopted a lower limit in the peak S/N ratio of 3 for the Gaussian fit to be accurate. Those locations in 
which the Gaussian fit proved difficult have been set to zero as for many positions in the [Nil] lines. As mentioned above, 
a necessary step is then the application of a correction for uneven slit illumination at every positi on, that varies between 
2 a nd 8kms~^ in our ca se. For a full discussion of this point and for details of the calculation see IBacciotti et al] (j2002[) 
and lWoitas et al.) (j2005D . Finally, the obtained velocity values were aggregated in a 2D map, which provides the velocity 
peak for each doublet and for each position resolved by STIS along and across both lobes of the jet. The radial velocity 
maps constructed in this way from the [O I] , [N II] , and [S II] doublets are shown in Fig [^. The maps present the same 
velocity asymmetry discussed in previous sections: the material in the blueshifted jet is in general twice as fast as in the 
redshifted lobe. The 2D maps also confirm that the [SII] (and also [Nil]) velocities at the base of the redshifted flow are 
lower than in [01]. Furthermore, the 2D maps allow us to confirm that in the approaching lobe the material with higher 
velocity tends to be located closer to the axis. Such peculiarities are discussed in Sect. 14.51 

Finally, we note that a marked asymmetry is present in the [O I] and [S II] images with respect to the position of 
the jet axis, especially in the red lobe. In fact, the material in the upper portions of the maps (i.e. right of the axis 
looking from the tip of the red lobe toward the source) has slightly more redshifted velocities than the material in the 
lowermost portion of the jet. This velocity asymmetry is known in stellar jets, and in previous papers has been tentatively 
interpreted as rotation of the jets around their axes (IBacciotti et"alll20'02t IWoitas et al.|[20b 5: Coff ev et al.ll2004 l2007l : 
IChrvsostomou et al.ll2008l ). The radial velocity maps of Fig. [5] are just another wa y of highlighting thi s effect for the 
RW Aur jet. The sense of the velocity asymmetry is consistent with that found in IWoitas et al. ll2005l) . and in other 



HST/STIS spectra of the same jet taken with the slit perpendicular to the jet axis (ICoffev et al.ll2004D . We add that, 
as mentioned in the Introduction, RW Aur is the only jet for which a discrepancy has been found between the rotation 
sense determined for the jet and for the accretion disk. The cause of this discrepancy has not yet been found, and a 
discussion of it is beyond the scope of the present paper. 



^ Colour version of the map is available in electronic form at 



http://www.edpsciences.org 
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3.4. Mass outflow rate 

In Paper I the average mass and momentum outflow rate Mj, Pj, respectively, have been estimated in an approximate 
way. Here we present a new, more detailed d e termi natio n of the Af; distrib ution along both jet lobes, exploiting two 
different methods described in iBacciotti et al.l (|1999[ ) and iNisini et al.l (|2005D . 

The first method combines the derived densities and radial velocities with estimates of the jet radius rj , to get the mass 
outflow rate through Mj = /xmH?T-H7rr?i;j. In this expression /j, = 1.41 the average atomic weight per H atom (lAllen fc Coxl 
120011) TTiH is the proton mass, and Vj is the bulk spatial velocity. Since tt-h is determined from line ratios, this method is 
not affected by reddening (as the wavelengths are so close) or distance errors. Since ratios formed with different elements 
are used, however, there is a dependence on the adopted elemental abundances. Furthermore, this method tends to 
overestimate the mass outflow rate because it assumes that the whole jet is uniformly filled at the derived total density, 
or, in other words, that the filling factor // is unity. This can be partially compensated for by the presence of regions in 
the flow even denser than those traced by the lines used in the BE technique. In any case, it should not be a significant 
factor affecting the RW Aur jet, as it appears to be very dense and coUimated. The bulk spatial velocities wj at each 
distance from the star were calculated by de-projecting th e radial velocities obtained previously with an inclination angle 
of the jet with respect to the line of sight of i = 46° ± 3 ()L6pez-Martfn et al.ll2003[ ). The radius of each jet section was 
assumed to be equal to one half the intrinsic F WHM {F WHMj^) of the jet brightness profile, determined from a jet 
image reconstructed from the seven spectra as in iPaper^ . after coadding all forbidden lines and velocity intervals, and 
correcting for the FWHM of the equivalent Point Spread Function (PSF) of the images, FWHMpsF- 

The uncorrected observed profile width, FVFTJMobs is completely analogous to the one obtained in iPaper ll for the 
first 2 arcseconds of the jet on both lobes, and it follows the same jet opening angle in the section from 2 to 4 arcseconds 
in the red lobe. Upon examining a reconstructed image of the stellar continuum, the FWHMpa,-p turned out to be 
O'.'l, i.e. comparable to the standard resolution limit of the instrument. The intrinsic jet width was then obtained as 
FWHM\^ = {FWHM^^^ — FWHAIpg^y^^ . The correction is larger for positions close to the star where the flow 
diameter is small, while it decreases with distance, becoming in practice negligible after 1'.'5 from the star. 

Combining these estimates with the total density tt-h reported in Figs. [5] and [3 obtained with the abundance set of 
lEsteban et all (|2004f ). one obtains A/j in both lobes at each position from the star, as illustrated in the bottom panels 
of the same figures, and shown in larger detail by the asterisks in Fig. [TOl In Figs. [5] and [71 one can see that the 
error bars of Mj are large in comparison with the range of this parameter. These errors are derived in an indirect way 

using Mj — iim}inYiTrr?Vj and combining errors for nn and wj into the mass flux errors according to the rules of error 
propagation. As a result, these error bars contain a systematic bias whereas the relative accuracy of the derived mass 
flux rates are much better and, in particular, the ratio of the rates in the red and blue lobes should be very reliable. 

The second method is based on the fact that the o bserved forbidden lin es are optically thin and thus their luminosity 
is proportional to the mass of the emitting gas (see iHartigan et al.l |1987[ ) . In STIS spectra the lines are calibrated in 
brightness units, thus one can convert them into absolute flux and thus luminosity units and estimate the mass outflow 
rate by comparing the total observed luminosity in one line with the value calculated theoretically. In this case Mj = 
jirriYi (tt-h V) Vt/h, with nn V = Lo\^s{^ine) / et\i{line), where V is the (unknown) volume effectively filled by the emitting 
gas, Vi is the tangential velocity of the emitting material, i.e. projected onto the plane of the sky. It is the length of the 
aperture considered, which we take equal to the extension of one pixel along the jet, Lobs(^*'T-e) is the observed luminosity 
in the li ne in the same regi on and et\i{line) is the emissivity per particle for the line calculated theoretically as explained 
in, e.g.. iPodio et al.l (|2006[ ). This method is affected by uncertainties in absolute calibrations, extinction, distance, and 
adopted elemental abundances, but does implicitly take into account the volume filling factor //, which, in practice, is 
also simply the ratio of Mj values derived from methods 2 and 1. For our calculation we used the total fiux of the [SII] 
and [O I] doublets and the brightest [N II] line at 6583 A. The adoption of the correct extinction value is quite critical for 
this method, as it can sensibly affect the observed luminosity and hence the mass loss rate directly. A local determination 
of Av with sufficient resolution is not available for this jet, and we have verified that this can introduce serious biases. The 
evidence, however, that the red lobe is, unusually, much brighter then the approaching lobe suggests that the extinction 
is different on the two sides of the system, and greater for the blue lobe, as if it propagate s into the the parent cloud. 
We have chosen therefore to apply to the two lobes the extremes of the range derived from lWhite fc Ghej (|200lh . that 
is, = 0.03 and 0.72 in the red and blue lobe, respectively. 

The values for the observed luminosity have then to be further corrected for the partial overlap of the slit area in 
the adjacent positions. The slit width is in fact O'.'l, but the positions are stepped by 0'.'07, which has the net effect of 
counting twice the intensity falling globally on an area of ^ O'.'IS times the slit length. We correct for this by multiplying 
the observed intensity by a factor 0.743 in both lobes. 

Finally, we take a distance of the system of 140 pc, in calculating ut, and et\^{line) is given as a by-product of the 
routines in our spectral diagnostic code. The mass outflow rates Mj calculated in this way are shown on a logarithmic 
scale in Fig. IIOI that illustrates the result determined from the various line luminosities and the comparison with the 
values obtained with the BE technique (flrst method, labelled 'from ratio'). The (quite signiflcant) uncertainties are 
similar for the two methods, and are reported in Figs. [5] and [71 and, as an average indicative value in Fig. [TO] with an 
error bar on the right hand side of the panels (A log Mq ^ 1.5 and 1.7 for the red and blue lobes, respectively). 

Overall, the mass outflow rate calculated with both methods is rather low, in the range 10~^-10~^ M© yr"-"^, and is 
almost constant along the beam. This is expected for a strongly coUimated outflow with no sideways dispersion as in this 
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case. Note however the decrease in Afj seen in the red lobe and corresponding to the faint region between knots J2 and 
J3. 

More precisely, excluding evident outliers, and, in both lobes, excluding the region beyond 2'.'1, that might be affected 
by stronger interaction with the environment, from the BE technique one derives on average Mj ~ 2.6 x lO^^MQyr"^ 
in the redshifted lobe and 2.0 x lO~^M0yr^^ in the blueshifted beam. 

The values estimated with the luminosity method show a similar distribution with distance from the star, but are 
higher by a factor 2-2.5 all along the jet, which indicates that probably an unidentified systematic effect comes into play. 
For example, the radius of the dynamically important material in the jet might be larger than 0.5 FWHM, as assumed for 
the 'ratio' method. Alternatively, the peculiar reconstruction of the jet images may introduce an unidentified geometrical 
effect in the calculation of the luminosity of the lines, possibly related to the unknown value of the local extinction. 
Because of the dependence of the latter method on extinction, distance of the object, and abundance/ionisation of the 
single species, we prefer to adopt the values derived from the BE technique. In any case the difference between the results 
from the two methods is much smaller than the uncertainty of the determinations. 

Overall, however, the results obtained regarding the mass outflow rates is very interesting, because it indicates that 
the jet has similar mass outflow rate on the two sides of the system, despite the evident asymmetries in velocity and 
density. This appears to suggest that the engine driving the jet is constrained on the two sides of the disk, as expected 
from axisymmetric magneto-hydrodynamic launch models. 

Further implications of this finding are discussed in the next Section. Here we just note that since approximately the 
same results are found with both methods, this suggests that indeed the filling factor is almost unity in this dense fiow. 
As mentioned above, however, this statement should be confirmed after a proper determination of the local extinction at 
each position, not available for the present study. 

4. Discussion 

In this paper we have derived with the BE diagnostic method the set of physical parameters characterising the emitting 
gas in both lobes of the RW Aur jet, with the best spatial resolution obtained so far along this fiow. These are illustrated 
in Figs. [5] to 1101 To ease the discussion and compare with previous estimates, we also give a summary of the obtained 
parameters in the form of average values within the first 2'.'1 along the jet in Tabled] 

In this Section we will discuss various aspects of the overall behaviour found in the derived quantities, starting from 
a comparison of the properties in the two lobes, in an effort to understand if differences in the ejection conditions may 
be present for the opposite sides. 

4.1. Comparison between the two lobes 

We consider the behaviour of the parameters in sections of equal length in the jet and counterjet, i.e. to distances of 
2'.'1 from the source. Note that this does not correspond to material ejected in the same elapsed time as we observe very 
different speeds on opposite sides of the central source. 

First of all, the electron density shows similar values in both the red and blue lobes, as well as the same gradual 
decrease with distance from the star, with oscillations corresponding to the variations of the total emission flux along the 
jet. The gradual decrease is consistent with a slow recombination of the g as from an initial ionisation state, combined 
with the gradual opening geometry of the jet, as discussed, e.g., in lBE99i The local increase of the electron density at 
the knot position can be justifled by the presenc e of shock fronts that generate locally a compression of the gas and an 
increase in the ionisation fraction (iHartigan fc Mors e 2002) ■ 

In contrast, the ionisation fraction Xc differs on average in the red and blue lobes. In the red lobe, Xc starts from 
very small values 0.02) and reaches its peak value at ^ 0.2. The ionisation fraction in the blue lobe at the positions 
closest to the source is higher 0.15) and reaches a maximum value of 0.4, which is two times higher than the peak 
value in the brighter red jet. The asymmetry of x^ in the two lobes has also been discovered in other bipolar jets (e.g. 
IIII30. Bacciotti et al.. 1999i) . where the fainter lobe similarly exhibits a larger excitation. The electron temperature is 
also higher in the blue lobe on average, (see Tabled]). This may indicate that the brightness of the jet lobes in this case 
depends at least in part on the excitation conditions within the flow, in addition to a peculiar extinction pattern (the 
approaching blue lobe appears to be less luminous). This may also be due to fact that the density is higher in the red 
lobe (see below). In both lobes one notes a gradual increase of Xc over the first arcsecond or so, and then a decreas e: such 
behaviour has been observed in several other cases at the beginning of jet fl ows (e.g. IIII30 and DG Tau. iBacciotti et al.l 
ll999tlLavallev-Fouauet et al.ll2000l : lBacciottiir2002l : lHartigan fc Mors^l2007t) . A related finding may be that in both lobes 
for positions close to the source the electron temperature To drops from 1.5 x lO^'K to about 10* K. This may give hints 
to the origin, still unknown, of the initial ionisation at the base of the jet. The behaviour of the temperature as a whole 
in the rest of the beam does not present any clear trend. It oscillates between 10* and 2 x 10* K, showing local increases 
at the position of the knots, but not always of the same amount. 

Finally, the derived hydrogen density tih in both lobes varies mostly in the range between 10* and 5 x 10^ cm^'^, and 
it seems to be higher in the red lobe than in the blue one - especially in the first arcsecond where it is higher by nearly 
one order of magnitude. This is in part the reason why the red lobe is brighter. 

It is worth noting that despite the higher density in the red lobe, the mass outfiow rate Mj is similar in both the red- 
and blueshifted jets. Note also that these rates were derived using the BE technique, which is independent of extinction 
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uncertainties. The fact that both lobes have approximately the same mass outflow rate, despite the numerous asymmetries 
in other quantities, may be important in modelling the launch mechanism (see below). 



4.2. Comparison with previous estimates 

Preliminary estimate s of som e quantities for the RW Aur jet observed with HST/STIS were presented in a previous 
work from our group (|Paper J ). A comparison of Fig. [1] of that paper with results presented here shows that the electron 
de nsities a re in good agreement, as well as the observed FWHM uncorrected for the FWHM of the equivalent PSF (Fig. 3 
m iPaper t ). even if our determination of this quantity comes from the coadded image in all forbidden lines instead of a 
single line. The mass outflow rates instead present discrepancies, which are discussed below. 

Tentative determination s of t he excitation parameters along the RW Aur j et hav e been attempted in two previous 
studies bv IPougados ct al] (j2002[ ) andLBacciotti et al. (1996). Doug ados et al.l (|2002[ ) quantified the physical properties 
from jet images obtained in 1997 at 0'.'2 angular resolution with the OASIS spectro-imager mounted at the Canada- 
France-Hawaii Telescope (CFHT). The authors derived the excitation conditions in the redshifted jet of RW Aur starting 
around 0'.'4 from the source using the same forbidd en lines (except [Nil] A6548), and diagnostic routines and elementa l 
abundances similar to the ones described bv lBE99l . The derived between 0'.'4 and 3" (Fig. 3 in lDougados et al.ll2002f) 
agrees with our results to within 20%. In contrast, the ionisation fraction Xq is generally smaller by about a factor of 
2 than the one calculated from our HST spectra, although similar behaviour with distance from the star is found in 
both papers. As a consequence the hydrogen density nu is also about twice the value determined here. Such differe nces 
in Xe niay come from the fact that in the present work a different abundan ce set has been assume d than in IBE99I (see 
the discussion in lPodio et al.|[2b06[ ). Nevertheless, the values found for To in iDougados et aTl ()2002f ) agree well with our 
estimates. Both papers report in the red lobe an increase in of up to 1.3-1.5x10^ K between 1" and I'.'S, followed by 
a drop down to 7000 K between I'.'S and 2". 

The excitation parameters in the red lobe were also analysed bv lBacciotti et all ()1996D . who used spectra with low 
spatial resolution (1'.'3) obtained with the Ca ssegrain twin sp e ctrogr aph of the 3.5-m telescope at Calar Alto Observatory 
in Spain i n 1993. The para. i neters derived by Bacciqtti et al.l (|1996[ ) present large discrepancies with respect to both our 



study and IDougados et al.l (j2002f ). iBacciotti et al.l ( 19961 ) found the ionisation fraction of the redshifted lobe decreases 



from a high value, Xe = 0.25, near the star to about 0.02 at 6-7". Within this distance a roughly constant hydrogen density 
of about 10^ cm~'^ and an extremely low temperature (typically 4500 K) were found. The reason for the big discrepancies 
is that i n that paper a version of the spectral diagnostic technique including the Ha line was used. As demonstrated 
later in lBE99l including the permitted Ha line can lead to incorrect determinations of the physical parameters, in 
particular higher ionisation fractions, and consistently smaller, unphysical temperatures. Contami nation by chromospheri c 
or magnetospheric Ha is also a cause of the high ionisation values found close to the source in IBacciotti etHI (|1996D . 
For these reasons the inclusion of Ha is now avoided in more recent diagnostic studies. The derived average mass outflow 
rate in IBacciotti et al.l (|l996l ). based on the results found at 4" from the source, was Mj ~ 5 x 10~®MQyr~-^. This is 
larger by one order o f magnitude than our values, basically because of the different jet radius adopted, i?j = 10^^ cm 
(IBacciotti et al.lll996l) . 

Recently, a set of HST/STIS spectra of jets from T Tauri stars including RW Aur, taken in 2002 close to the so urce 
with the slit perpendicular to the flow have been analysed with the BE diagnostic technique bv lCoffev et al.l (j2008l ). Of 
the two spectra, taken at O'.'S from the source, only the one for the red lobe presented enough signal in the [Nil] lines to 
be diagnosed. This ga .ve a Position (with respect to the jet axis) - Velocity diagram in the physical quantities, illustrated 
in Fig. 15 and 16 of ICoffev et al.l (|2008( ). The average electron density is 2 x 10"'cm~'^, the ionisation fraction varies 
between 0.02 and 0.07, and the total density is around 3 x lO^cm"'^. These values are consistent with those found in 
the present work a t the same position from the star. The electron temperature, however, turns out to be only 5 x lO'' K 
(jCoffev et al.ll2008l ) a factor of 2 lower than observed here. This might be due to a very low value of the [S II]/[0 1] ratio, 
which introduces large uncertainties in the temperature determination. In the same paper, an estimate has been given 
of the mass outflow rate in the flow, that turns out to be 1. 7 x 10~^ Mp, yT~^, twice as large as our estimate at the same 
position. The difference comes mainly from the fact that in lCoffev et al.l (|2008l ) the FWHM measured in [01] was larger 
by a factor of 1.4 than measured by us because no correction for the FWHM of the PSF was applied. Such a correction 
is necessary particulary close to the star (see Sect. 13. 4| ). 

The mass outflow rate determined in the presen t study f or the r ed lobe is al so smaller by a factor 2-5 (depending on 
position along the jet) than the one determined in iPaper 1 and in iHartigan et al. (199H) from the spatially unresolved 
[O I] line luminosit y at the source. This d iscrepancy is due to the fact that in [Paper ll the ionisation fraction was taken 
from the study of IPougados et aTl (|2002| ). which, as mentioned above, reports smaller ionisation fractions, and hence 
higher total densities than o urs. In a ddition, in the absence of any other determination, the same ionisation fraction was 
adopted for the blue lobe in lPapei^ . while we found here that the average Xc in the approaching lobe is more than twice 
that in the red lobe, which reduces the derived total densities. In addition, no correction for the PSF width was applied 
in iPaper 1 which explains their larger values of Mj in both lobes close to the source. 

The reported differences may also be due to time variability of the ejection properties of the source, as the epoch of 
the observations analysed in the various studies differ (the data were taken in 1993, 1997, 2000 and 2002). T Tauri stars 
are known to vary rapidly in luminosity, and new jet knots have been seen to arise from jet sources on timescales of one 
year or less. 
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Finally, it is interesting to estimate from our measurements the ratio of mass ejection to mass accretion rates Mj/Afacc, 
and compare it with those for MHD mo dels for jet launch. An accretion rate of A/acc = 1.6 x 10~^ Mq yr~^ was determined 
for RW Aur by iHartigan et al.l ()1995D from measurements of the veiling of photospheric lines caused by the accretion 
shock. For the ejection rate, one should consider the sum of the mass outflow rates of both jet beams. Adopting our 
average mass flux rate one obtains Mj = 4.6 x 10~^ M(7) yr~^, th en i\/j/Macc ~ 0.003, which is compatible with MHD 
models with large values of the magnetic lever arm (|Ferreira et a l. 2006). Note however that jValenti et al. (1993) also 
analysed the veiling in RW Aur spectra and derived a lower accretion rate of A'/acc = 3.4 x 10^^ Mq yr^^ which suggests 
that Mj/Macc may be 5 times higher. Non-stationary accretion could also be at play, as RW Aur is k nown to have 
a highly variable spectrum , in which emission lines undergo strong and rapid changes fe.g. iMundt fc Giampapail98^ 
IStout-Batalha et al.ll200(3l ). This has been ascribed to fast variations of the accretion rate. 



4.3. Comparison with other young stellar Jets 

The same set of physical parameters (including the total hydrogen density) have been derived with O'.'l spatial resolution 
for the DG Tau and HH30 jets, as well as for the jet from the Herbig Ae/Be star LkHa233 from HST images and 
sp ectra. This allows us a direct comparison with the present results for the RW Aur jet. According to the findings 
in iBacciotti et ahl (|1999[ ). the total hydrogen density of HH30 averaged over the first A" of i ts redshifted jet is about 
3 x lO^cm-3, while for LkHQ;233 it is 1.1 x 10^ cm'^ within the first 2" (jMelnikov et al.ll2008D . Thus, these jets appear 
to have lower densities than the RW Aur jet, for which th e average is lO^cm"'^. In contrast, the DG Tau jet reaches 
values as high as 10^ cm~'^ at its base (IBacciotti et al.|[2002l ). similar to the first positions in the RW Aur red lobe. Other 
parameters in these j e ts present similarities with RW Au r. As in the RW Aur red lobe, in both HH30 and DG Tau 
(jBacciotti et al.l l2000l 120021 : iLavallev-Fouauet et al.l l2000t ) the ionisation fraction is seen to increase in the initial jet 
section, reach a plateau at 1-2 arcsec from the source, and then decrease slowly following a recombination curve. This 
has been interpreted as the ionisa tion fraction being produced close to the source, and then being 'frozen' into the jet 
medium. IHartigan fc Mor"s3 (|2007t ). however, demonstrated through their analysis of 2-epoch STIS slitless spectra that in 
HH 30 there are ionisation peaks that follow the movement of the knots. Thus at least part of the ionisation is produced 
and maintained locally at the internal shock fronts. This is consistent with what we observe in RW Aur, where we find 
a correlation between the peaks of rto, Xc and To and the positi ons of the bright features. This is particularly evident at 
knots 32, J3, and J6 in the red lobe. Furthermore, as noted bv IHartigan fc Morsel (j2007f ) in HH30, we find a tendency 
of the ionisation fraction to increase upstream of the knot, inste ad of a t the knot itself. A similar effect is also seen in 
the region before the bright condensation A in the HH 46/47 jet (|BE99t ). 

Several HH je ts have been investigated on larger scales using ground-based telescopes (|BE99l : iPodio et al.l 120061 : 
iNisini et al.ll2005[ ). The hydrogen density of knots in these jets is generally found to be between 600 and 6 x lO^cm"^, 
that is one to three orders of magnitude lower than what we measure in the first arcscconds of the RW Aur jet. This 
indicates not only that the RW Aur jet is one of the densest measured so far, but also that the material in the jet tends 
to be diluted as it propagates to large scales, both by beam opening and because of sideways losses due to internal bow 
shocks. 



4.4. Proper motion of the knots 

The distribution of forbidden line emission along the RW Aur jet reveals several knots clearl y visib le in the 'total 
flux-d istance' diagram. Positions of these knots have also been reported in iLopez-Martm et al.l (j2003D and iPvo et all 
(j2006t ). who investigated the propagation of some bright knots along the jet. With this goal in mind, they c ompared 
ground-based obser v ations taken in 1997-2002, with our HST data taken in Dec. 2000 and first published in iPaperl . 
iLopez-Martm et al.l (j2003[ ) found that the mean proper motion of knots for the red- and blueshi fted jets are 0'.'16 and 
0'.'26yr~^, respectively, and this agrees with the velocity asymmetry of the flows. IPvo et al.l (|2006l ) traced also the knots 
in [Fell] Al. 644 /xm. Comparison of knot positions in the red jet in the IR [Fell] line with those in optical forbidden 
emission de rived from HST spe c tra sh ows that the knot at 0'.'25 exists at the same position in all investigated data 
sets and so. iLopez-Martm et al.l (|2003D and lPvo et al.l (I2006D sugge s t that the knot may be stationary. In the blue lobe. 



there is another knot visible in the [Fell] at 0'.'17 which Pvo et ahl (|2006( ) suggest is also stationary. At the same time. 



the HST data do not reveal any maxima at those positions in the optical. Instead, the forbidden emission flux steadily 
decreases from the position closest to the central source. The lack of the emission maximum in the optical might be due 
to problems with the extraction of the continuum as the star is very bright in this spectral range. 

Using the pr oper motion of O '.'l 6 yr^^ proposed for the red jet (jPvo et al.ll2006[ ) we suggest that the [Fell] maxima 
at 1'.'2 and 2'.'29 (jPvo et al.ll2006l) correspond to knot J2 (0'.'9) and J4 (1'.'95), respectively, allowing for the 2-year period 
between the observations. While the positions of [Fell] knots in the blue jet at 0'.'48 and 1'.'25 coincide well with knots 
All (0'.'5) and A12 (1'.'2) from the HST data, this may be coincidental. Instead taking into account the hig her proper 
motion in the blue lobe, the knot All at 0'.'5 (which is quite wide) may correspond to the [Fell] knot at 1"25. iPvo et al.l 
(l2006l) came to the same conclusion that the shoulder around 1'.'25 may correspond to All at 0'.'5 (= knot B4 from 
iLopez-Martm et al.ll2003l ) if the knot has a proper motion of 0'.'34±0.1 yr~^. 
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4.5. Jet asymmetry 

One of the most peculiar aspects of the RW Aur jet is that it shows as ymmetries in sever al respects. First of all, the jet 
lobes show considerable difference in their velocities, noted already in lHirth et al.l (jl994l ). The average peak velocity of 
the lines in the blue lobe is about — 190kms~^, while a velocity of only +100kms~^ is measured on the oppo site side. 
This asymmetry persists close to the source, as confirmed by the first analysis of our STfS dataset in iPaper ll . Here we 
can add that the [01] and [SII] lines exhibit a systemic shift of velocities visible at the beginning of the red lobe: the 
velocity in [SII] is 10-15 km s~^ lower than in [01]. The fluctuations of radial velocity along each lobe, however, are less 
than 15-20%, with the largest variation seen in the prominent knots J3-J4. In the red beam the velocities scatter around 
+100kms~^, whereas in the blue beam a very slight decrease can be noted. 

Additional evidence of the velocity asymmetry comes from analysis of the for bidden line s in ot her spectral ranges. 
Radial velocities of the [Fell] A1.66/^m line of RW Aur jet have been reported bv iPvo et al.l (|2006t ). Unlike the optical 
forbidden lines, the radial velocities of the infrared [Fell] in the redshifted jet steadily increase from about lOOkms"^ 
at the beginning of the lobe to about 140kms~^ at 3". Because [Fell] traces denser regions of the jet than [SII], this 
implies that the denser material may have had a higher velocity in the past, or is accelerating, whereas the velocity of 
material traced by optical [SII] Unes remains almost constant. This suggests that denser material is located closer to the 
axis, where higher velocities are predicted by theory. On the other hand this effect is not seen in [O I], which should also 
trace higher densities. In the blue jet the [Fe II] velocities behave similarly to the optical forbidden lines: they present a 
gradual decrease with distance from the jet source in the same range of velocities (—200 170kms~^). 

The ionisation fraction x^. also shows evident asymmetries, while the electron density and temperature are similar in 
both jet lobes (see Table [T]). In the blueshifted lobe the average Xc is 0.23, while in the redshifted one it is 0.08. As a 
result the redshifted jet is found to be 2 times denser. On the other hand the redshifted jet is about 2 times slower and 
this makes the mass outflow rates similar on the two sides. 

To flnd an explanation for the observed asymmetries is an open problem. The fact that the mass outflow rate is about 
the same in both lobes suggests that the engine driving the jet has certain constraining symmetries on both sides of the 
disk. Therefore, the origin of the differences in density, velocity, and excitation may instead reside in an inhomogeneity 
of the ambient medium around the star or, in a magneto-centrifugal launching scenario, in a different configuration of 
the circumstellar magnetic field that shapes both the accretion funnels and the open field surfaces along which the jet is 
launched. 

The question if such an asymmetric configuration can be generated and remain stable for a sufficiently long time needs 
to be addressed with detailed theoretical modelling. It is known that a number of young stars do show asymmetries in the 
magnetic field distribution on their surfaces. The magneticall y generated cold spots that cover a co nsiderable portion of the 
surface of weak-line T Tauri stars (WTTS) (see, for example. lGrankiiilllQQStlGrankin et al.ll2008f ) pre sent Doppler images 
consi stent with largely inhomogeneous distributions and different properties in opposite hemispheres (jRice fc Strassmeieij 
[19961). Stron g asy mmetries in the dis tribution of spots hav e also been revealed in evolved late-types stars, e.g., XX Tri 
( Strassmeieij|l999t ). At the same time. lGrankin et al.l (|2008[ ) from analysis of the periodic variability of WTTS, argue that 
in some cases exte nded spotted regions can be stable over a timescale of at least two decades (see the WTTS V410 Tau in 
IStrassmeieij|1999l ). Recently, a complex non-dipolar topology of the stellar magnetic field has been discovered in t he tw o 
accreting T Tauri stars V2129 Oph and BP Tau, through spectro-polarimetric observations ijDonati et al.l [20071 . |2008[ ). 
However, a connection between surface magnetic properties and large-scale circumstellar magnetic configuration seems 
unestablishcd for the time being. On the other hand, th ere are theoretical models that propose a large-scale asymmetrical 
topology for the circumstellar magnetic field (see, e.g.. lLong et al.|[2007T ) which may ultimately explain the asymmetry 
of jet properties. 

5. Conclusions 

In this paper we analyse in detail the first few arcseconds of both lobes of the RW Aur jet, using long-slit high-resolution 
HST/STIS spectr a. The physic al properties of the outflow are derived applying the so-called BE technique to the spectra, 
described, e.g., in lBE99l and in lPodio et al.l (|2006f ). Our results can be summarised as follows: 

1. The distribution of forbidden emission along both outflows reveals several local maxima identified as jet knots. We see 
seven emission knots in the redshifted lobe within 3'.'8 from the source, and at least two knots in the blueshifted lobe 
within 2'.'1. The distribution of emission within the knots in the blue lobe appear to be flatter than in the redshifted 
section. 

2. The electron density ric in both outflows is close to the high density limit (log(no)ciit = 4.39) in the region near the 
driving source and decreases gradually with distance from the star. The electron temperature in the redshifted 
lobe decreases steeply within the first 0'.'3 from 1.5 x 10"* K to 10^ K and then varies around this value. The same 
trend is observed in the first few positions in the blueshifted lobe, but then Tc reaches about lO^'K and subsequently 
shows higher values with a maximum of about 2 x 10^ K. 

3. We found that the gas in the red lobe has an ionisation fraction Xc ranging between 0.02 and 0.2, with an average 
of 0.08; in the blue lobe this parameter is higher and varies between 0.1 and 0.4 (average 0.23). x^ increases in the 
first section of both lobes, reaches a maximum at 1-2 arcseconds from the source and then starts to decrease. This 
behaviour appears to be a common feature also for other j ets from young stars. Our results for Xe for the red lobe 
agree within a factor of two with those of previous studies (jPougados et al.ll2002l : [Bacciotti et al.lll996t ). 
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4. The hydrogen density rin in the brighter redshifted lobe appears to be twice as high (mean value nn ~ 6.5 x lO^cm"'^) 
than in the fainter blueshifted section, w here an average value of uh ^ 2.3 x 10* cm"'^ is found. RW Aur jet is found 
to be the second densest (after DG Tau, iLavallev-Fouauet et al.|[2000D of the T Tauri jets investigated so far. 

5. The kinematical properties of the flow have been fully re-analysed at high resolution and illustrated in accurate 
graphic detail. We confirm that the asymmetry in radial velocity (— 190kms~"'^ and -|-100kms~^ in the blue and red 
lobes, respectively) persists in all emission lines and down t o the source, and t hat the flow presents indications of 
rotation around its axis, in the same sense as determined in lWoitas et al.l (|2005f ). The pattern of proper motions of 
the knots has been updated with the data presented here. 

6. Strong variations of ne,Xe and Te are observed close to knots J3-J4 in the red lobe and A12 in the blue lobe. The 
radial velocity variations are very limited and reach a maximum of 20% of the average at knots J3-J4 in the red lobe. 
Variations of the physical parameters at the other knot positions are suggested but are much less prominent. 

7. The average mass flux Mj calculated from the forbidden line ratios is similar in the first 2'.'1 of both lobes: 2.6 x 
lO~^M0yr~^ for the red- and 2.0 x 10~^MQyr~^ for the blueshifted jet. The Mj/Macc ratio is rather low and we 
estimate it to be between 0.003 and 0.012, depending on the adopted values for the mass accretion rate. This ratio, 
however, can still be accommodated by magneto- hydrodynamic models having a large value of the lever arm. We note 
however that the ratio may also be influenced by non-stationary accretion. 

8. The overall trend of the physical parameters along the first few arcseconds of the RW Aur jet is similar to that of 
HH 30 and DG Tau (|Bacciotti et al.lll999l l2000l l2002l : iHartigan fc Morsel l2007l) . which have been investigated with 
similar resolution. Analogies in the properties of the RW Aur, HH 30, DG Tau jets in regions close to central star 
can refiect analogies in the mechanisms operating in that region, suggesting the same engine is accelerating the jets 
in the T Tauri stars with outfiows. 

9. Our study of the RW Aur jet indicates for the first time that, despite the detected marked asymmetries in physical 
and kinematic properties between the two lobes, the mass outflow rates in the two lobes are similar. This appears 
to indicate that the central engine has constraining symmetries on both sides of the system, and that the observed 
asymmetries are probably due to different environmental conditions. We discuss the likelihood that these are related 
to a large-scale inhomogencity in the distribution of the local magnetic fields. 

In conclusion, our HST/STIS study of the physical parameters of the bipolar jet from RW Aur has highlighted a 
number of properties that appear to be common to many objects of the same class when observed close to their source. 
That said, it reveals that peculiar elements of this and other jets, like the observed asymmetries in the two opposite 
lobes, remain unexplained even after analysis conducted at sub-arcsecond resolution. Nevertheless, this work shows that 
an accurate investigation close to the source of dynamically relevant properties, like the mass outflow rate, that present 
no asymmetry, can constrain models. Obviously whether mass outflows rate is balanced in other asymmetrical bipolar jets 
should be investigated. The combined efforts of future theoretical studies, and soon-to-become available interferometric 
observations at milliarcsecond resolution are expected to give useful signposts to an understanding of these phenomena. 
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Table 1. Average physical parameters within the first 2'.'1 along the RW Aur bipolar jet lobes obtained with the 
application of the BE technique. 



Parameters redshifted blueshifted 

Tie (cm--') 3600 4100 

([SII],[OI]A6363) (kms-^) 100 -180 

0.08 0.23 

Tc (K) 12 300 16 400 

nn (cm-^) 65 900 22 900 

Mj (Mq yr-i) 2.6 x 10"® 2.0 x 10"" 
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RWAur jet [SII] X6731 




-2 2 -2 2 -2 2 -2 2 -2 2 -2 2 -2 2 

Vj. [100 km s'] 



Fig. 1. Position- Velocity diagrams of the [S II] A6731 line, obtained for the seven positions of the STIS slit stepped across 
the RW Aur jet. The slit was oriented parallel to the jet with P.A. = 130 deg. All velocities are systemic in this and all 
subsequent figures (with v^, = +23.5kms~^ in the heliocentric frame). The contours are drawn logarithmically within 
the same range of surface brightness in all the panels (floor: —14.7, i.e. 2 x 10~^^crgs~^ cm~^ arcsec"^ contour step: 
0.2 in logarithm). 
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RW Aur jet 

[OI] 6300 [OI] 6363 [Nil] 6548 [Nil] 6583 [SII] 6716 [SII] 6731 




[100 km s"^] 



Fig. 2. Position- Velocity diagrams in each of the forbidden lines in the RW Aur jet, obtained from coadding the spectra 
in the seven slit positions across the jet. High spatial resolution (O'.'l) is retained in the direction along the jet. The 
logarithmic contours are drawn from —14.2 (6.3 x 10~^^ergs~^ cm~^ arcscc"^ with incremental logarithmic step of 

0.2. 
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RW Aur redshifted lobe 
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Fig. 3. Ratios of the optical forbidden lines along the redshifted lobe of the bipolar jet from RW Aur. From top 
to bottom: intensity profile resulting from the sum of forbidden lines integrated in velocity and across the jet width, 



in units of 10 ergs 



^, and with labelled knots; [S II] A6716/A6731 ratio and its error, with superim- 



posed theoretical values in the low (dotted line) and high (dot-dashed line) density limits, respectively, calculated 
for Te = 10^ K; [O I] A6300/A6363 with its error, compared with the theoretical ratio of Einstein's A coefficients 
(dotted line); same as above, for the [Nil] A6583/A6548 ratio; [Nil] A(6548+6583)/[OI] A(6300-|-6363), with its error; 
[01] A(6300-f6363)/[SII] A(6716-f6731), with its error. The positions where corrupted values were corrected (see text) 
are marked by filled symbols. 
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Fig. 4. Same as Fig. [S] but for the blueshifted lobe of the RW Aur jet. 
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Fig. 5. Calculated values of the integrated flux in [OIJA6300 (asterisks) computed as Flux([0 1] A6300)caic = 
Flux([OI] A6363)obs x 3 according to the theoretical ratio of Einstein's A coefficients, are plotted together with the ob- 
served integrated fluxes of [01] A6300 (line). This demonstrates that no significant flux is lost because of the [01] A6300 
line lying at the left border of the detector. 
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Fig. 6. Physical properties of the emitting gas in the redshifted lobe of the RW Aur jet as a function of distance from the 
source, derived from line ratios in Figs. [3] and 3] with the BE technique. From top to bottom: intensity profile as in Fig. [31 
radial velocity in the [01] A6300 and [SII] A6731 lines, integrated across the jet (uncertainty between 3 and 30kms~^, 
see text); electron density Uc] hydrogen ionisation fraction Xc', electron temperature Te; total density tih = n^/xc and 
mass outflow rate Mj calculated from the above total density with the BE technique (see Section 13.41 and asterisks in 
Fig. [101). 
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Fig. 7. Same as Fig. |51 but for the blueshifted lobe of the RW Aur jet. 
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Fig. 8. Radial velocity v, measured in the [01] AA6300,6363, [SII] AA6716,6731, and [Nil] AA6548,6583 lines emitted by 
the RW Aur jet, as a function of distance from the central source. Top panels: distribution in the first (2'.'1) of the two 
lobes. Bottom panel: w,. for the whole red lobe section (3'.'9) analysed. [O I] A6300 is marked by empty squares, [O I] A6363 
- asterisks, [SII] A6716 - triangles, [SII] A6731 - diamonds, [Nil] A6548 ~ small filled squares and [Nil] A6583 ~ pluses. 
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Fig. 9. Spatial distribution of peak absolute radial velocity \vr\ in the lobes of the RW Aur jet, derived by aggregating 
the values obtained by Gaussian fits to the line profiles in each of the seven spectra and at varying distance from the star. 
From top to bottom: \vr\-Taa.ps obtained from the sum of the [SII] lines; same for the [01] doublet (only the [01] A6363 
line was considered in the blueshifted lobe because of the proximity of the [01] A6300 line to the CCD left edge); same 
as above, for the [Nil] doublet. Positions where the Gaussian fit failed (because of low S/N) have their value set to zero. 
The colour scale for |?Ji-| ranges from to 210kms~^. All velocities are corrected for the heliocentric radial velocity of 
RW Aur A (+23.5kms-i). 
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Fig. 10. Comparison of mass outflow rates Mj in the two lobes of the RW Aur jet calculated from the total density 
obtained with the BE technique (method I, asterisks, labelled 'from ratio') and derived using absolute fluxes of forbidden 
lines (method II, diamonds: [S II] A(6716+6731), triangles: [O I] A(6300+6363), and squares: [NII]A6583. The typical 
uncertainty on Mq for each lobe is also showed with an average error bar at the right hand side of the panels. The 
averages correspond to A log A'/© ~ 1.5 for the red- and 1.7 for the blueshifted lobe. 



